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Abstract 

By employing the Bianchi identities for the Riemann tensor in 
conjunction with the Einstein equations, we construct a first order 
symmetric hyperbolic system for the evolution part of the Cauchy 
problem of general relativity. In this system, the metric evolves at zero 
speed with respect to observers at rest in a foliation of spacetime by 
spacelike hypersurfaces while the curvature and connection propagate 
at the speed of light. The system has no unphysical characteristics, 
and matter sources can be included. 
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1 Introduction 



We consider the evolution part of the Cauchy problem in General Relativity ||] 
as the time history of the two fundamental forms of a spacelike hypersurface: 
its metric g and its extrinsic curvature K. On such a hypersurface, for ex- 
ample an "initial" one, these two quadratic forms must satisfy four initial 
value or constraint equations. These constraints can be posed and solved as 
an elliptic system by known methods that will not be discussed here. (See, 
for example, [[!]].) 

The Ricci tensor of the spacetime metric can be displayed in a straightfor- 
ward 3+1 decomposition giving the time derivatives of g and K in terms of 
the space derivatives of these quantities. These expressions contain also the 
lapse and shift functions characterizing the threading of the spacelike hyper- 
surfaces by time lines. However, proof of the existence of a causal evolution 
in local Sobolev spaces into an Einsteinian spacetime does not result directly 
from these equations, which do not form a hyperbolic system for arbitrary 
lapse and shift, despite the fact that their characteristics are only the light 
cone and the normal to the time slices 0. 

In this paper we present in detail, using "algebraic gauge" or a gener- 
alized harmonic time slicing condition on the lapse function, a hyperbolic 
system for these geometric unknowns that is based directly on the Riemann 
tensor, using the Bel|§ electric and magnetic fields and inductions, and the 
full Bianchi identities. Its characteristics are the physical light cone and the 
time direction orthogonal to the spacelike submanifolds M t of the chosen slic- 
ing. The spacetime curvature tensor is among the quantities that are evolved 
explicitly. Its propagation is coupled to that of the connection, extrinsic cur- 
vature, and metric of M t through relations introduced by H. Friedrich0], who 
obtained an analogous vacuum system based on the Weyl tensor, causal but 
with additional characteristics. Our system constitutes a complete first or- 
der symmetric hyperbolic system equivalent to the usual Einstein equations 
including matter. This new system, which was sketched briefly in H , is com- 
pletely equivalent to our previous one || || 0. The space coordinates and the 
shift three- vector are arbitrary; and, in this sense, the system is independent 
of a gauge choice. In fact, we can also make the lapse function take arbitrary 
values by introducing a given arbitrary function into its definition. Indeed, 
our procedure for the lapse function is equivalent to specifying g~ x ^ 2 N as an 
arbitrary positive function of the spacetime coordinates. 
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2 Metric, Connection, and Curvature 

It is convenient to write the metric on the hyperbolic spacetime manifold 
M t x R as 

with 9° = dt and 9 l = dx l + (3 l dt, where (3 l is the shift vector. The cobasis 
6 a satisfies 

de a = —cr^Aff 1 , (2) 

with C l oj = — C % jQ = dj/3 % and all other structure coefficients zero. The 
corresponding vector basis is given by e = d t — /3 3 dj, where d t = d/dt, 
di = d/dx\ and the action of eo on space scalars is the Pfaffian or convective 
derivative 

e [f] = d f = d t f-^d 3 f. 

We shall assume throughout that the lapse function N > and the space 
metric g on M t is properly Riemannian. Note that g^ = gij and g lJ = g 13 . 
(An overbar denotes a spatial tensor or operator.) 

The connection one-forms are defined by 

^ Q /3 7 = r a ,3 7 + g a5 C e s(/3g y )e - ( 3 ) 

where (^7) = \(f3^ + 7/5), [^7] = \{f3^ — 7/9), and T denotes a Christoffel 
symbol. We adhere to the convention 

VaV P = dy + V°^ ca (4) 

for the spacetime covariant derivative V Q and likewise for the spatial covari- 
ant derivative Vj. The Riemann curvature tensor is given by 

R a p,pa = d p iU a f3 a — 8 a UJ a fj p + u" \ p UJ X p a — id" \ a UJ X /3 p — UJ a ' (3\C X pa (5) 

and the corresponding Ricci identity is 

V a V/3M 7 - V/3V Q M 7 = R a(3 ^ S U 5 . (6) 

(The comma in R a p,<y5 is used to distinguish the two antisymmetric index 
pairs [afl] and [7 5].) 



3 



The connection coefficients are written in 3+1 form as 





y\% T"^ 

J- jk J- j'/c; 












= NgVdjN, 


U Oi 


= u° l0 = N- 1 d l N, 


^ 00 


= N^doN, 



(7) 



and 
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, l -^N- 2 d g ij = -N- 1 K ij , (8) 

where is the extrinsic curvature (second fundamental tensor) of M t and 
for any t-dependent space tensor T, we define another such tensor <9oT of 
the same type by setting, as in 



4 = ^ - C 0} (9) 

where Cp is the Lie derivative on M t with respect to the spatial vector (3. 
Note that Bq and di = d/dx l commute. 

3 Bianchi Equations 

We recall that the Riemann tensor of a pseudo-Riemannian metric satisfies 
the Bianchi identities 

V a i2/3 7i A/L( + V 7 -R Q/ 3 i AAt + V ' pR^a^n = 0. (10) 

These identities imply by contraction and use of the symmetries of the Rie- 
mann tensor 

V a R a pM + V M i? A/3 - V x Rtf = 0, (11) 
where the Ricci tensor is defined by 

If the Ricci tensor satisfies the Einstein equations 

Ra/3 — Pa/3, (12) 
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then the previous identities imply the equations 



^7aR a p,x^ = V A p M /j - (13) 

Equation (|TDJ) with {af3^} = {ijk} and ( |T3| ) with (3 = do not contain 
derivatives of the Riemann tensor transversal to M t ; hence, we consider these 
equations as constraints ("Bianchi constraints"). We shall consider the re- 
maining equations in fllUD and ( |13"D as applying to a double two-form A a/3 ^, 



which is simply a spacetime tensor antisymmetric in its first and last pairs 
of indices. We do not suppose a priori a symmetry between the two pairs 
of antisymmetric indices. These equations, called from here on "Bianchi 
equations," read as follows 

VoA hkM + V k A QhM + V h A k0M = 0, (14) 

V A\am + V h A h iM = V x p^ - V M p A * = J\/iit (15) 

where the pair [X/j] is either [Oj] or [jl]. We next introduce H two "electric" 
and two "magnetic" space tensors associated with the double two-form A, 
in analogy to the electric and magnetic vectors associated with the electro- 
magnetic two-form F. That is, we define the "electric" tensors by 

En = A° ifij = -N~ 2 A 0lflj , (16) 

Dij = —TjihkVjlmA ' , 

while the "magnetic" tensors are given by 

Hij = ^N^mA^m, (17) 



Bji = 2^ lr 1ihkA j hk . 



In these formulae, rjijk is the volume form of the space metric g. 

Lemma. (1) If the double two-form A is symmetric with respect to its two 
pairs of antisymmetric indices, then = Ej i; = Dj i; and = B^. 
(2) If A is a symmetric double two-form such that 

A a = A x aj xf3 = cg af3 , (18) 
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then Hij = Hji = Bji = Bij and Eij = . 

Proof. (1) If A is a symmetric double two- form, the proof is immediate. 
(2) The Lanczos identity [§] for a symmetric double two-form, with a tilde 
representing the spacetime double dual, is given by 

Aap^ + A a/3: x^ = C a \gp^ — Ca^gpx + Cj3^g a \ — Cp\g afl , (19) 

with 

C af3 = A a(3 -^Ag af3 , (20) 
A = A Xa M . 

The equalities E = D and B = H follow by a simple calculation that employs 
the relation ^oijfc = Nrjijk between the spacetime and space volume forms. 

In order to extend the treatment to the non-vacuum case and to avoid 
introducing unphysical characteristics in the solution of the Bianchi equa- 
tions, we will keep as independent unknowns the four tensors E, D, B, and 
H, which will not be regarded necessarily as symmetric. The symmetries 
will be imposed eventually on the initial data and shown to be conserved by 
evolution. 

We now express the Bianchi equations in terms of the time dependent 
space tensors E, H, D, and B. We use the following relations found by 
inverting the equations ( f[6l) and (|17|) 



Aoi, 0j = ~iV%, (21) 

A hk , oj = NrfhkHij, (22) 

Ahk,lm = rfhkrf ImDij, (23) 

A 0j ,hk = NrfhkBji. (24) 



We will express spacetime covariant derivatives V of the spacetime tensor A 
in terms of space covariant derivatives V and time derivatives do of E, H, 
D, and B by using the connection coefficients in 3+1 form as given in (0) 
and ©. 

The first Bianchi equation (|i~4D with [Xfj] = [Oj] has the form 

^oAhk,oj + ^kAoh,oj — ^hAokflj = 0. (25) 
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A calculation incorporating ([?]), (|) and (^1])-(^4|), then grouping derivatives 
using do and V«, yields the first Bianchi equation in the form 

doiv'hkHij) + 2NV [h E k]j + (LJmj = 0, (26) 

{L x ) hk j = NK l jV i hk H il + 2{V [h N)E k]j + 2Nr ] i lj K l [k B h]i (27) 
-(V^rA^'A- 
The second Bianchi equation (|15D, with [A/i] = [Oj], has the form 

Vo^4°j,oj + Vfti j^j = Jojj) (28) 

where J is zero in vacuum. A calculation similar to the one above yields for 
the second Bianchi equation 

doEij - N^iVhHtj + (L 2 )y = J oji , (29) 

(L 2 ) tJ = -AT(trK)^ + NK k jE ik + 2NK l k E k] (30) 
-{VhN)^ 1 ^ + NK k hV lh iV m kj D lm + (V k N) V l kj B u . 

The non-principal terms in the first two Bianchi equations fl2"E| ) and ([19]) 
are linear in E, D, B, and H, as well as in the other geometrical elements N~K 
and VN. The characteristic matrix of the principal terms is symmetrizable. 
The unknowns Ei(j\ and H{u-\, with fixed j and % = 1, 2, 3 appear only in the 
equations with given j. The other unknowns appear in non-principal terms. 
The characteristic matrix is composed of three blocks around the diagonal, 
each corresponding to one given j. 

The j th block of the characteristic matrix in an orthonormal frame for 
the space metric g, with unknowns listed horizontally and equations listed 
vertically, (j is suppressed) is given by 







Et 


Ei 


E 3 




H 2 


H 3 




(29)x 


/ 


£o 











m 3 




\ 









Co 





-m 3 






























(P%3 







-N£ 3 




£o 










(26) 3 i 














Co 







(26) 12 
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-N£ 2 
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/ 



(31) 
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This matrix is symmetric and its determinant is the characteristic polynomial 
of the E, H system. It is given by 



-N b (^e)(u a r- (32) 

The characteristic matrix is symmetric in an orthonormal space frame and 
the timelike direction defined by 8q has a coefficient matrix T that is positive 
definite (here T is the unit matrix). Therefore, the first order system is 
symmetrizable hyperbolic. We do not have to compute the symmetrized 
form explicitly because we will obtain the energy estimate directly by using 
the contravariant associates E % \ H % \ ... of the unknowns. 

The second pair of Bianchi equations is obtained from ([HP and (|T5| ) with 



[A/i] = [lm]. We obtain from (|i~4 ) 



doWhkV 3 imDij) + 2Nr] J i m V [k B h]j + (L 3 ) hkj i m = 0, (33) 

[L 3 )hkM = 2Nr ] n j[m K j l] r) i hk D in + 2r) j lm (V [k N)B h]j (34) 
+2NK l[h E k]m + 2NK m[k E h]l + 2H i[l (V m] N)r t i hk . 



Analogously, from uHh we obtain 

WimBij) - N v kh ^ lm V h D kj + (L 4 ) lM = -NJ lmi , (35) 

(L 4 ) tM = -N(tiK)r l h m B i3 + 2Nr l h 3[m K3 l] B ih + 2NK h rf lm B h] (36) 
-(V j N)n hi in n lm D hn - 2Nrf hi H j[m K h l] + 2E i[m V l] N. 



Consider the system ( ^3|) and (|35D with [lm] fixed. Then j in rjji m is also 
fixed. The characteristic matrix for the [lm] equations, with unknowns Di 3 
and Bij, j fixed, with an orthonormal space frame, is the same as the matrix 

©• ' 

If the spacetime metric g is considered as given, as well as the sources, 
the Bianchi equations (fffip, (P^D, (PH[), and ( P3| ) form a linear symmetric 
hyperbolic system with domain of dependence determined by the light cone 
of g. The coefficients of the terms of order zero are ViV or iVK. The system 
is homogeneous in vacuum (zero sources). 
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4 Bel Energy in a Strip 

Multiply (H) by \r)i hk H l i and recall that r)i hk rf hh = 25\, r]irkV ihk = ^i5 h r - 
5\5 h u and d g ij = 2NK ij . Then we find that 

^HVdoWhkHij) = ^doiHijlfS) - M 1; (37) 

M x = ^iK.H ln rf hk H i jd Q (^g ka 9 im ) (38) 
= iV((trK)iP J ' - if^^ + Ki j H il )Hij. 
Likewise, multiply (p5| ) by to obtain 

^'4^- = ^o(^^') - M 2 , (39) 

M 2 ee N(K\E lj + Ki j E il )Eij. (40) 

Multiply Q33D by (l/4)r} r h % lm D rs and (^) by (l/2)rj h lm B ih to obtain anal- 
ogous results for the second pair of Bianchi equations. Sum the expressions 
so obtained from the four Bianchi equations (|26l), (p9|) , (|33|) , and (|35|) . The 
spatial derivatives add to form an exact spatial divergence, just as for all 
symmetric systems. Indeed, we obtain 

^ (|E| 2 + |H| 2 + |D| 2 + |B| 2 ) + V h {NE i 'rf h i H lj ) (41) 

-VhiNB^^Du) = Q(E, H, D, B) + S, 

where we have denoted by | • | the pointwise g norm of a space tensor, and 
where Q is a quadratic form with coefficients ViV and iVK given by 

Q = -'l ni h ''H"(L 1 ) kKj - £«(L 2 ) (j - - A nr hk ri,' m D"{L !s ) hk , m (42) 

+M 1 + M 2 + M 3 + M 4 . 
The source term S, zero in vacuum, is 

SEE^-iiVW"^. (43) 



9 



We define the Bel energy at time t of the field (E,H,D,B), called a 
"Bianchi field" when it satisfies the Bianchi equations, to be the integral 

B(t) = U (|E| 2 + |H| 2 + |D| 2 + |B| 2 ) Mgi . (44) 

Z J Mt 

We prove the following 

Theorem 1. Suppose that g is C 1 on M x [0,T] and that the g t norms of 
ViV and NK. are uniformly bounded on M t; t e [0, T]. Denote by n(t) the 
supremum 

7r(0=Sup Mt (|VJV| + |JVK|). (45) 

Suppose the matter source J G L 1 ([0, T], L 2 (M t )) ; then the Bel energy of a 
C 1 Bianchi field with compact support in space satisfies for < t < T the 
following inequality 

B(t) < |s(0) + ^||5|| L2{ M T )rfr|exp(C J\(r)dr), (46) 

where C is a given positive number. 

Proof. We integrate the identity ( f4"T|) on the strip M x [0,t] with respect 
to the volume element Hg T dr. If the Bianchi field has support compact in 
space, the integral of the space divergence term vanishes. The integration 
of a function d$F on a strip of spacetime with respect to the volume form 
dr fMg T goes as follows 



f f doFdr^ =11 {d t -(3 l d i )Fdr HT 

JO JM T ' JO JM T 

= [* f d T (F(detg T ) 1 / 2 )drd 3 x (47) 

JO J Mr 

- f Q J m (F(det g T )-^ 2 d T (det g.) 1 / 2 ) /ig T dr 

- f f ViifiF) fi gr dr+ f [ FV^ Ht dr. 

Jo Jm t ' Jo J M T 

Therefore, if F has compact support in space, we can express the right hand 
side of (0) as 

/ F Hr -[ Ffa+fl N{tiK)F^ T dr, (48) 
J M T Jm Jo Jm t 
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where we have used 



d T (detg T )V 2 = -(detg^A^trK) + (detg^V^ (49) 



with trK = K 3 j. 

The integration of (41) on a strip leads therefore to the equality 

B(t) = B(0) + f / (Q + S) /i gT dr, (50) 

J JM T 

with 

Q = Q + iiV(trK)(|E| 2 + |H| 2 + |D| 2 + |B| 2 ). 

We deduce from this equality and the expression ( 42|) for Q the following 
inequality, with C some number 

Bit) < B(0) + Cj* Sup Mr (|ViV| + \NK\)B(t) dr + J* \\S\\ L 2 (Mt) dr. (51) 

This inequality and the resolution of the corresponding equality imply the 
result. 

Corollary. If the metrics gt are complete, the inequality given in the theorem 
extends to C 1 Bianchi fields that are square integrable in /xg t . 

Proof. One uses a truncating sequence. 

Remark 1. The quantities — VfciV and —2K{j are, respectively, the (Ok) 
and (ij) components of the Lie derivative £ n g of the spacetime metric g 
with respect to the unit normal n of M t . [Its (00) component is identically 
zero.] The Bel energy, therefore, is conserved if this Lie derivative is zero. 

Remark 2. Energy estimates for the double two- form A can be deduced 
from the Bel-Robinson tensor of A when A is assumed to be symmetric in its 
pairs of indices and when a timelike vector field has been chosen0. However, 
we do not make such a hypothesis about the symmetry of A here. 



5 Local Energy Estimate 

We take as a domain Q of spacetime the closure of a connected open set whose 
boundary dfl consists of three parts: a domain u t of M t , a domain ou of M , 
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and a lateral boundary L. We assume L is spacelike or null and "ingoing," 
that is, timelike lines entering at a point of L are past directed. We also 
assume that the boundary dQ is regular in the sense of Stoke's formula. We 
use the identity (flip previously found and integrate this identity on Q with 
respect to the volume form \x^ T dr. Let f(t, x) = t+(p(x) be the local equation 
of L. Set Vi = dif = di4> and z/ = 1 — /JV*; these are the components in the 
coframe 8 a of the spacetime gradient d a f = v a . Then we have from (fTT|) 



\f (|E| 2 +|H| 2 +|D| 2 + |B| 2 )/ig t -i / (|E| 2 + |H| 2 + |D| 2 + |B| 2 )/i g0 



' ' — 12 i |tt|2 i |T-v|2 i |T3|2\ , at-,, (Vci a U\h /r> A T\\h 



+ 



L 



-z/ (|E| 2 + |H| 2 + |D| 2 + |B| 2 ) + Nu h {(E A H) ft - (B A D) A } 



f/ {Q + S) Hr dT. (52) 



We have set 



(E A H) h = ^i^Htj = £(E A H)£ 



00' 



that is, (EAH)^, for each fixed j, is the vector product in three dimensions 

of the vectors and Hu-\. Therefore, the g norm of the vector (E A H) h 
satisfies 

|E A H| < \E U) A H U} \ < £ \E {j) \\H U) \ < ^(|E| 2 + |H| 2 ). (53) 

Let v = (z/j), vq > 0, and \u\g < N~ x vq, that is, let v_ = (uq, Pj) be timelike 
or null, as it must be for L. [Note that our sign conventions for y_ are suitable 
for applications of Stokes's theorem in a spacetime with signature ( — h++).] 
We can now deduce from flSBp and an analogous inequality for (D, B) that the 
integral over L in (1521) is non-negative. It is strictly positive if u is timelike. 



Remark. The integral over a null boundary L is also strictly positive if ([53]) 
is a strict inequality. This will be the case if any of the norms of the vector 
products of E^) by jj^ j s \ ess than the product of the norms of these vectors, 
that is, if any of these pairs of vectors are non-orthogonal: E^H^ ^ or 
D l ^B i(j) ^ for some j. An analogous property has been obtained through 
the use of the Bel-Robinson tensor of the Weyl curvature in ||T0 |. 
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We define B(u>t) by replacing in the definition of B(t) the integral over 
M t by an integral over uu t 

B(uh) = U (|E| 2 + |H| 2 + |D| 2 + |B|> gt . (54) 

Z JuJt 

A result of the non-negativity of the integral on L is that the inequality ([51]) 
holds on Q with £>(r) replaced by B(u T ), uj t = M T fl Q, and M T replaced by 
uj t . Therefore the inequality (|46|) also holds with the same replacements. In 
particular, we have B(uj t ) = if B{ojq) = and J = (vacuum case). Then 
E = H = D = B = Oinfiif they vanish on the intersection of Mo with the 
past of Q. (Cf. for related results |10| and Note that such a result is 



not sufficient to prove the propagation of gravitation with the speed of light 
because it treats only curvature tensors that are zero in some domain, not 
the difference of non-zero curvature tensors. The Bianchi equations are not 
by themselves sufficient to estimate such differences because their coefficients 
depend on the metric, which itself depends on the curvature. 

In the next section, we will give a further first order symmetric hyperbolic 
system linking the metric and the connection to our Bianchi field. This 
further system is inspired by an analogous one constructed in conjunction 
with the Weyl tensor by Friedrich ||. 

6 Determination of (f , K) from Knowledge of 
the Bianchi Fields 

We will need the 3 + 1 decomposition of the Riemann tensor, which is found 
by combining (^), (0), and fl|); namely, 

Rij,ki = Rij,ki + 2Ki[kKi]ji (55) 
Rajk = 2NV u K k]i , (56) 
Roiflj = NidoKij + NKikK^ + VifyN). (57) 
From these formulae one obtains the following ones for the Ricci curvature 

Ra/3 — R 1 a,if} 

Rij = R ij -N- 1 d K ij + K ij txK-2K ik K k j -N- 1 Vid j N, (58) 

R Qj = NidjtvK-VhK^), (59) 

R 00 = N(d trK- NKijK ij + V%N). (60) 
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The identity 

Bogy = -2NK l3 (61) 
and the expression for the spatial Christoffel symbols give 

4r\ = V h (NKij) - 2V (i (iVK i) fc ). (62) 

Therefore, from the identity fl56|) , we obtain the identity 

dot % + NV h K, 3 = K l0 d h N - 2K h {l d j} N - 2R 0(hj) h . (63) 

On the other hand, the identities fl5"7| ) and (E3) imply the identity 

doKij + NRij + VjdiN ee -2NR° ifij - iV(trK)i^ + NRij. (64) 

We obtain equations relating F and K to a double two-form A and matter 
sources by replacing, in the identities ( |63"|) and (|64]), Roa,j) by (Aorjj) + 
^(i,«)o)/2, -R°i,oj by (^4°i,oi + ^ j,oi)/2, and the Ricci tensor of spacetime by 
a given tensor p, zero in vacuum. The terms involving A are then replaced 
by Bianchi fields as in (0)-(|24]). 

The first set of identities ( |63"D leads to equations with principal terms 

d t h ij + Ng hk d k K ij . (65) 



To deduce from the second identity fl64l ) equations which will form together 
with the previous ones a symmetric hyperbolic system, we set 

AT = «- 1 (detg) 1/2 , (66) 

where a is a given scalar density of weight one. This is the "algebraic gauge" 
[ |T2]j . (It can also be considered as a change in the name of the unknown from 
N to a [0.) The condition (|66|), if t denotes the Christoffel symbols of g, 
implies that 

t h ih = d i ]ogN + d i ]oga. (67) 



The second set of identities ( p4[ ) now yields the following equations, where 
iV denotes a _1 (det g) 1 / 2 , 

doKij + Ndht^ = N^iht^m-it^-diloga^jk-djloga)] (68) 
+N(d i d j log a - t k ijd k log a) - N(E l3 + E jt ) - iV(trK)^ + N p iy 
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The first set (|63| ) yields 

BoT^j + NV h Kij = NK ijg hk (T m mk -d k loga) (69) 
-2NK h (t (T m j)m - dj) log a) - Nr} k { *B j)k - NH^rfJ. 

We see from the principal parts of (|B"8| ) and ( |69"D that the system obtained 
for K and F has a characteristic matrix composed of six blocks around the 
diagonal, each block a four-by-four matrix that is symmetrizable hyperbolic 
with characteristic polynomial N 4 (£o£°)(£ a £ a ). The characteristic matrix in 
a spatial orthonormal frame has blocks of the form 



( ^ 






N£ 3 


\ 




£o 















Co 







\ N£ 3 








6> 


/ 



(70) 



Remark. In the equations considered above, N is to be replaced by a _1 (det g) 1 / 2 . 
Instead of this algebraic replacement, we can consider an evolution equation 
for N of the type 

d Q N + NhrK = Nf, (71) 

where / is an arbitrary function on spacetime. For / = 0, (|71|) is the harmonic 
slicing condition. For arbitrary /, the general solution of ([H]), if g and K 



iY = 7 - 1 (detg) 1 / 2 , (72) 



are linked as in (61), is 
with 7 such that 

4 7 + fry = 0, (73) 

that is, 7 satisfies a linear first order differential equation depending only on 
the known quantities (3 and /. It does not depend on any of the previously 
defined unknowns (E, H, D, B, g, K, F) of our system in algebraic gauge. The 



modified system consists of the equations in algebraic gauge [(|26|), (p9|), (|33|). 
(|35D, (|6l|), (|68|) , (|69D 1 but with iV now replaced by ([72]) . Therefore, a is 
replaced by 7. The additional equation is (|7"T|), or, on account of (|72|), this 
role can be regarded as played by (ff3|). 
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7 Symmetric Hyperbolic System for u = (E, H, 
D,B,g,K,f; 



We denote by S the system composed of the equations (|26|), (p9[), (|33|), (|35|), 
(|61~1), (|68|), and (|69|) , where the lapse function iV is replaced by ct _1 (det g) 1 / 2 . 
This system is satisfied by solutions of the Einstein equations whose shift (3, 
hidden in the operator do, has the given arbitrary values and whose lapse 
has the form iV = a _1 (det g) 1//2 . (Clearly, any N > can be written in this 
form.) The following lemma results from the previous paragraphs. 

Lemma. For arbitrary a and (3, and given matter sources p, the system S 
is a first order symmetric hyperbolic system for the unknowns u. 

Note that in this Lemma, the various elements E, H, D, B, g, K, and 
f are considered as independent. For example, a priori, we neither know 
that F denotes the Christoffel symbols of g nor that E, H, D, and B are 
identified with components of the Riemann tensor of spacetime. 

An analogous Lemma holds for the system S' that has an additional 
unknown N ( equivalent ly, 7) and an additional equation ([n]) [equivalently, 



(|73|)1 with / replacing a as an arbitrary function. 
We now consider the vacuum case. 

Initial Values. The original Cauchy data for the Einstein equations are, 
with (ft a properly Riemannian metric and tft a second rank tensor on Mq, 

g|o = eft, K| = ift- (74) 

The tensors (ft and ift must satisfy the constraints, which read in vacuum, 

Roi = 0, (75) 

= G o = R 00 -\gooR (76) 

= # 00 + ijVV / X 3 

= IiV 2 (i?-^.^ + (trK) 2 ), 

with R = g^Rij. The initial data given by (ft determine the Cauchy data 
T h ij\ and thus -R^wlo- Then, Rij^i\o is determined by using also ift. To 
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determine the initial values of the other components of the unknown u of 
the system S, we use the arbitrarily given data a and (3. In particular, we 
use N = a _1 (det g) 1 / 2 to find Roi,jk\o an d to compute VjdiN\ appearing in 
the identity (|58|). We deduce from ( 58[) d Kij\ when = 0, which enables 
Roi.oj to be found from fl57|). All of the components of the Riemann tensor of 
spacetime are then known on Mq. We identify them with the corresponding 
components of the double two-form A on Mq: the latter have thus initially 
the same symmetries as the Riemann tensor. We find the initial values of 
(E, H, D, B) according to their definitions in terms of A. 



8 Existence Theorems 

In order to define Sobolev spaces on a manifold M, it is convenient to endow 
M with a C°° Riemannian metric e. We will suppose that e has a non-zero 
injectivity radius, 5q > 0, hence is complete, and that it has Riemannian 
curvature uniformly bounded on M, as well as its derivatives of all orders 
relevant for the required Sobolev properties to hold. 

With I an interval of R, and under the hypotheses made on e, each 
manifold Mxl endowed with the metric e—dt 2 is globally hyperbolic |15|, [13 



We denote by D the covariant derivative associated with e (the "e-covariant 
derivative"). 

Theorem 2. Hypotheses, a. The arbitrary quantities are such that, with I 
an interval containing 0, a is continuous on M x I , and there exist numbers 
Ai, A 2 > such that on M x I 

< At < a < A 2 , 

Ba,{3eC°(I,H s+1 ). 
b. The initial data are such that 

D0, ^ € H s+1 , 

and <p is a continuous properly Riemannian metric on M uniformly equiv- 
alent to e. That is, there exist strictly positive numbers b\, b 2 such that at 
each point of M and for any tangent vector £ to M at this point, 
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Conclusion. If s > 3, there exists a number T > 0, [0, T] G I, such that the 
system S has one and only one C 1 solution u on [0, T] x M taking the Cauchy 
data u deduced from <p, if), a, (3. The solution u is such that the components 
ofu different from g, as well asg — (f), are in C°([0, T], H s ) nC 1 ([0, T], H 8 -i). 
There exist B x and B 2 such that 

3ig(*,*)<g(6£) £)■ 



Proof. 1. We replace the equation ((6T|) in the system S by the equation 

4 ^ 3 - <f>ij) = ~2NK i:j + Cpfay (77) 

We still have a symmetric first order system S for the unknown u, deduced 
from u by replacing g by g — <fi. This quasilinear system is hyperbolic as 
long as g is properly Riemannian. 

The hypotheses made on the arbitrary data a, (3 and on the initial data 
d), ifi imply, through the Sobolev multiplication properties, that the initial 
data u|q belong to H s whenever s > 1. 

2. If the Riemannian manifold (M, e) is the Euclidean space R 3 , the existence 
theorem is known, with s > 3/2 + 1. [Linear case: Friedrichs ||16|| . Quasilinear 
case: Fischer-Marsden [pi , using semigroup methods, and Majda (T^|, using 
energy estimates.] 

A proof under the hypotheses we have made on (M, e) and on the data 
can be obtained along the same lines, using energy estimates, because the 
spaces H s on (M, e) have the same functional properties that they have on 
R 3 . The energy inequalities for u and its spatially covariant e-derivatives 
of order < s can be obtained directly on [0, T] x M, as we have indicated 
in the case of the Bianchi equations. One uses s-energy estimates for linear 
symmetric first order hyperbolic systems and the Garding duality method 
l!J to prove existence and uniqueness for the linearized system. One uses 



the contraction mapping principle in the norm C ([0, T], H a -i) on the one 
hand, and boundedness in the norm C°([0, T], H s ) fl C 1 ([0, T], H s _i) on the 
other hand, as in Majda @, to prove the theorem for our quasilinear system. 
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Remark. The spacetime defined by the g component of u together with a 
and (3 on M x [0,T] is globally hyperbolic, with a C 1 metric if a and /3 are 
C 1 . 

Domain of dependence. The exterior sheet of the characteristic cone of 
the system S at a point of spacetime V is the light cone of the spacetime 
metric g. It has been proven by Leray |T5|], for globally hyperbolic Leray 



systems on a manifold, that the values of a solution of the Cauchy problem 
in a compact set uj t C M t depend only on the values of the Cauchy data in 
the compact set lj that is the intersection of M with the past of uj t - We 
will prove this result for our symmetric hyperbolic system. 

We recall that on a spacetime (V, g), V — M x R, temporally oriented 
by the orientation of R, the past V(x) of a point x is the union of the paths 
timelike with respect to g and ending at x. The past V(u>) of a subset u C V 
is the union of the past of its points. 

Lemma. Let (V, g) be a globally hyperbolic manifold with V = M x R and 
g a C 1 metric. Let ujt be a compact subset of Mt with non-empty interior 
ujt and boundary Ouit- Then, 

q = v{u T ) n {t > 0} 

is a compact subset of V. Its boundary is 

dVl = uo T u uj U L. 

The lateral boundary L is generated by null geodesic arcs issuing from points 
ofdu>T- At an interior point of such an arc, L has a null tangent hyperplane 
and is ingoing. 

Proof. On a globally hyperbolic manifold, the past of a closed set is closed; 
and the intersection of the past of a compact set with a "past compact set" 
(here, t > 0) is compact. Therefore, fl is compact. 

If y & Q H {0 < t < T}, there exists x G ojt such that y G V(x). Suppose 
that the path xy is strictly timelike. Then there is an open neighborhood 
V(y) of y in V(x) such that z G V(y) implies that z G V{x). Consequently, 
y cannot be on the boundary L, which must therefore be generated by null 
paths. Analogous reasoning shows that x cannot be an interior point of ujt- 
For the proof of the rest of this Lemma, see Leray fl5 l. Lemma 97. 
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Recognizing that in a C 1 globally hyperbolic manifold the lateral bound- 
ary L admits a null tangent hyperplane almost everywhere, we prove the 
following theorem. 

Theorem 3. (Domain of Dependence) Two C 1 solutions u x and u 2 of S 
coincide on ujt if the arbitrary data a, (3 coincide in the past of ujt defined 
by gi, while the Cauchy data (</> 1 ,i/? 1 ) and ((fi 2 ,'if) 2 ) coincide on u , which is 
the intersection of the gi-past of lot with Mq. 

Proof. The difference of two solutions satisfies a linear homogeneous first 
order symmetric hyperbolic system with C 1 principal coefficients and C° non- 
principal ones. We set 7 = l\ — r 2 , k = Ki — K 2 , G = gi — g 2 , a = Ai — A 2 , 
and likewise for the Bianchi fields obtained from A. The linear system is 
composed of three sets of equations in which we denote by I various linear 
forms in 7, k, . . . with coefficients analytic in \\\ and u 2 . 
The first set is 

d^n + N^duh, = 4 7) , (78) 
doh. + NA^^ = if. (79) 
The second set consists of an equation for e = E x — E 2 

4ey - N^g^ninJm = 4 e) (80) 

and similar equations for h = Hi — H 2 , d = Di — D 2 , and b = Bi — B 2 
deduced from the Bianchi equations. [Note that e in (|8CiD is not to be confused 
with the metric e introduced in Sec. 8 and below in the discussion of Hf A \] 
The third set of equations is 

4<?ii = 4 G) - ( 81 ) 

Raising and lowering indices with the metric gi, we multiply ( [78| ) by 7^, 
( |79"D by (|80P by e lJ and so on. We add the results and integrate over a 
domain Q of the type indicated in the preceding Lemma, constructed with 
gi. We denote by 8{uJt) the integral 

E{u t ) = \( (| 7 | 2 + |k| 2 + |e| 2 + |h| 2 + |d| 2 + |b| 2 + |G| 2 )/i gl . (82) 

Integration of (^) over Q leads to the identity 

£{ut) = £{u ) + [ cr + [ [ q/j,g (t) dt, (83) 

JL JO Ju>t 
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where cr is a three-form and q is a quadratic function of 7, k, . . . , G with 
coefficients analytic in ui and 112. The three-form cr arises from integration of 
a d derivative together with a spatial divergence from the principal terms of 
our equations. We made an analogous computation for the Bianchi equations 
and showed that the integral over L in that case is positive. The same 
argument holds here for the terms coming from the second set of equations. 
It is easy to prove that the integral over L is also positive for the first and 
third sets of equations. On the other hand, q is bounded on Q, up to a 
multiplicative constant C, by the sum of the squares of the gi norms of 
7, k, e, . . . , G. Hence, we have the inequality 



It follows that £{ijJ T ) = if S(u ) = 0, which is the stated result. 

By using Theorem 3, we will show that the local existence result holds for 
a non-compact manifold M without any hypotheses on the fall-off at infinity 
of the initial data. 

Definition. A tensor u on M belongs to a space Hf h ("u.l." denotes 
"uniformly local") if there exists a covering C of M by open sets uj^ with 
compact closures such that 

1. The restriction ofu to each of these open sets belongs to H s {uj^), that 
is, u is square-integrable on u)u\ in the volume element defined by the 
metric e introduced in Sec. 8, and so are all its e-covariant derivatives 
of order < s. 



Note that different choices of uniformly locally finite coverings, as well as 
different uniformly equivalent metrics e, define the same Hf L space because 
they define equivalent norms. 

Lemma. If, for the covering C, the arbitrary quantities and the original 
initial data are such that 




T 



(84) 



2. 



Su Pill u ll// s (^ w ) = \\ u \\h^- < 00. 



a > 0, 
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f3eC°(I,H^ 1 ), 

and defines a properly Riemannian metric on M , uniformly equivalent to 
e, then the initial data for u belong to H^ L if s > 1 . 

Proof. The proof follows from the definition and the Sobolev multiplication 
properties. 

Theorem 4. Hypotheses, a. The arbitrary quantities are such that, with I 
an interval containing 0, 

aeC°(I,H^ 2 ), a>A>0, 

/3e£ (/,# s +i)- 
b. The initial data are such that 

and 4> defines a properly Riemannian metric on M , uniformly equivalent to 
e. 

Conclusion. If s > 3, there exists a number T > 0, [0,T] C J, such that the 
system S has one and only one solution u G C°([0, T], Hf ■) nC 1 ([0, T], Hf_\) 
taking the Cauchy data u G i?"' 1 ' deduced from 0, i/j. 

Proof. Suppose M is not compact. Consider a covering of Mt by compact 
sets iiTj . . The pasts of these sets in the globally hyperbolic metric —A 2 dt 2 + 
B 2 e. cover Vr = M x [0,T]. (The numbers A and B 2 here refer to Theorem 
2.) The intersections of these pasts with M define a covering of M by 
compact sets Kq\ There exist truncating functions that are smooth, 
bounded uniformly with respect to i ( cLS £1X6 their e-covariant derivatives of 
order < s), equal to one on K$ and equal to zero outside another compact 
set strictly containing Kq^ Use the sets to define the H£ l 

norms of tensors on M. 

Consider the Cauchy problem for the system S with Cauchy data Uq = 
^''Uo- There exists a constant C independent of i such that 

|| Uq \\h s < C||Uo||i?u.i.. 
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Corresponding inequalities hold for the given quantities a, (3. Theorem 2 
shows that there exists a number T' < T, independent of i, such that S has 
one and only one solution on M x [0, T] that takes the initial value Uq . 

Remark. The associated metric gW is uniformly equivalent to e and the 
spacetime metric gW satisfies the inequalities posited in Theorem 2. There- 
fore, the past of a subset in this metric is contained in the past defined with 
e. 

We obtain, from the set of iz^'s defined on Vt, a solution u for the Cauchy 
problem for S on V t by setting 

u(x,t) = u®(x,t) on V(K^) n {0 < t < T'} 

because, by Theorem 3 and the previous Remark, we have 

u (j \x,t) = u®(x,t) on V{K$) n V(K$) n {0 < t < T'}. 

Finally, we prove 

Theorem 5. The solution of S with the data deduced from the arbitrary 
quantities a, f3 and the initial data <f>, if) is a solution of the vacuum Einstein 
equations with lapse N = a _1 (det g) 1//2 and shift (3 if the initial data satisfy 
the constraints. 

Proof. The result holds for H s+ \ as well as for H%+\ data. We spell it out for 
the first case. It has been proven by Choquet-Bruhat and York||, for initial 
data in the indicated function spaces and satisfying the constraints, and for 
a given arbitrary shift (3, that the vacuum Einstein equations have one and 
only one solution in the "algebraic gauge," that is, with a lapse function in 
the given form, such that g 6 C°(Mx [0, T]), (Dg, K) G C°([0, T],H S+1 ). The 
connection of g, the second fundamental form K and the Riemann tensor of g 
satisfy the system S and take the Cauchy data Uo- They coincide, therefore, 
with the solution u found previously. In particular, the component g of u 
together with TV = a _1 (det g) 1 / 2 and (3 are a metric that satisfies the vacuum 
Einstein equations. 

Remark. It is a result of the coincidence of the components g and K of u 
with the solution of the Einstein equations with initial data (f> and if), given 
shift (3, and lapse iV = a~ 1 (det g) 1 / 2 that these components are in fact such 
that (Dg, K) G C°([0, T],H s+ x). ' 
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